The exponential parameterization of the quark mixing matrix 
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We comment on the exponential parameterization of the quark mixing matrix, by stressing that 
it naturally incorporates the Cabibbo structure and the hierarchical features of the Wolfenstein 
form. We extend our results to the neutrino mixing and introduce an exponential generator of the 
tribimaximal matrix. 
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I. INTRODUCTION 

The quark mixing matrix can be written in different 
ways, any of the proposed forms displays nice features 
and disadvantages. Whatever form one uses, four arbi- 
trary parameters and the assumption of its unitarity are 
necessary to get physically meaningful results. The mod- 
els can be roughly grouped in two categories, the first 
inspired to Euler like rotation matrices, the second, con- 
taining explicit hierarchical features, employs an expan- 
sion, around the unit matrix, in term of some key parame- 
ters. The original Kobayashi and Maskawa matrix[l| had 
been written in terms of three mixing angles #1,2,3 and 
one CP violating phase 8. In this parameterization the 
first family decouples from the others in the limit 9± — > 0. 
The particle data group chooses a form in which the 
CP violating term is appended to the matrix entries re- 
sponsible for the coupling of the first and third genera- 
tions of quark mass eigenstates. Finally Wolfenstein [3| 
has proposed a matrix emerging from a kind of pertur- 
bative expansion in terms of the Cabibbo coupling pa- 
rameter A = 0.22[4|. A third model, bridging between 
the ([l|,0|) and [J is based on the so called exponen- 
tial parameterization, which emerges from the request of 
unitarity, automatically satisfied by setting Q 



V =e A 
At = —A 



(1) 



The second condition in Eq. ([I]), expressing the anti- 
hermiticity of the matrix, is ensured by the following 
specific choice 



except A3, are real. In the spirit of Wolfeinstein crite- 
ria, we use the Cabibbo strength A as key parameter and 
make the following identifications [6| 



Ai =A 
A 2 =y\ 2 
A3 =x\ 3 e iS 



(3) 



containing an implicit hierarchical assumption on the 
coupling between the different quark families. The van- 
ishing of the x, y coefficients allows the decoupling from 
the b sector reducing the matrix to the s-d Cabibbo mix- 
ing, namely 



V, 



(z,S/)->0 



/ 



V 



0\ 


0) 



cos(A) sin(A) 
— sin(A) cos(A) 
1 



( 4 ) 

It is also to be stressed that V^->o = lj which means 
that the parameterization in ([3]) contains the assumption 
that the vanishing of the Cabibbo parameter determines 
the decoupling of the entire quark matrix. The phase 8 
is associated, as in the particle data group choice, with 
the smallest coupling term. In this paper we will see 
how the exponential parameterization yields a flexible 
tool to analyse the quark mixing phenomenology and the 
relevant consequences. 



A : 



Ax A 3 
— Ai A 2 

-AJ5 -A 2 



(2) 



The vanishing of the diagonal entries secures that the 
matrix V be unimodulaiQ- The sub-labels 1, 2, 3 deter- 
mine the mixing d-s, s-b, d-b respectively, all the entries, 
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1 It would be sufficient to have a matrix with null trace, but for 
practical reasons we use the form(2). 



II. THE MATRIX A AND THE WOLFENSTEIN 
PARAMETERIZATION 



We will prove that the quark mixing matrix written as 
in Eq. ([1]) naturally contains the Wolfenstein parameter- 
ization and the Euler like forms as well. By keeping the 
expansion of the exponential in Eq. ([1} up to third order 
in A, namely 



V=t+A+ 



A 2 



A3 
3! 



A^ 
4! 



o(A 5 ) 



(5) 



2 



we obtain the mixing matrix in the form 
1 - 



2 



A 
4! 



F S 



^4GA 3 



A 

A* , A 
ACA 2 



I _ Al + Al _ 
1 2+4! 



(AA 



AFX 3 

ABA 2 



1 - 



A A 



A = y, F = p — it], G = 1 — p — irj 
B = 2 — A 2 (p - i - ttj), B + C = -2\ 2 {p - i) 
P = fcos(5) + i, »/ = -2Bin(<y) 



(6) 

Eq. © is recognized as a Wolfenstein-type parameteri- 
zation, the Taylor expansion at higher order can provide 
more accurate expansion in the Cabibbo coupling param- 
eter, as we will see in the following. The expansion at the 
third order allows a one to one correspondence between 
the Wolfcnstcin parameters and those of the matrix A, 
which can be written in the form 







A = 



-AX 3 {p 



IT] 



h) 



A AX 3 (p — irj 
AX 2 
-AX 2 



(7) 

Using for A, p, r\ the following values, close to those 
given in the literature Q: 



A = 0.2272 ±0.0010, 



a - 221 + 064 
p — u.zzi_ 028 , 



A = 0.8181°;°?? 
V = 0.34018:82 



we hnd for x and 6 the following values 



-0.359l8;°i, 6 = 0.88318;^ 



and we get for the mixing matris0 
\A\ = 



0.97429 0.22523 3.86 • 10" 3 
0.22512 0.97341 0.04215 
8.10- 10~ 3 0.04154 0.99910 



(8) 



in good agreement with the values reported in [7]. 
Higher orders expansions will be considered in the forth- 
coming sections. 



III. THE GEOMETRICAL MEANING OF THE 
EXPONENTIAL PARAMETERIZATION AND 
THE EULER LIKE FORMS 

We have so far proved that the exponential parameter- 
ization of the mixing matrix has some nice features which 



This result has been obtained by expanding the matrix at any 

arbitrary order, namely V = and by keeping N=50. We 

o n - 

have not included the errors deriving from the experimental and 
systematic uncertainties, the relevant analysis would require ex- 
treme care for fitting the data and such an effort is out of the 
purposes of the present note. 



makes its use quite interesting. Before going further let 
us speculate on the geometrical (physical) meaning of the 
matrix A, which can be understood as a kind of Hamil- 
tonian ruling the process of quark mixing. We introduce, 
therefore, the Schroedinger equation 

id T f_ = %> (9) 
where i))\ t —q are the quark mass eigenstates, and 

H cx lA (10) 

Within such a picture the matrix V is the evolution op- 
erator associated with Eq. ©. In the case of vanishing 
CP phase 6^0, the Hamiltonian in (fl"0| can be written 
in terms of SO (3) generators, namely 



H = XRi+ yX 2 R 2 + xX 3 R 3 



(11) 



with 





1 
















• 


, R2 = i\ 





■ 









\o 


-1 





Ri 



1 
R 3 = i | 

-10 



The Schroedinger equation §§§ can, accordingly, be 
viewed as a vector equation of the type 



it x 



(12) 



A(-yA,2;A 2 ,-l) 



where = (tpi, ^2^3) is the vector associated with the 
quark held. The problem of the quark mixing is therefore 
understood as a rotation, induced by an Euler-like torque 
equation. 

The torque vector £1 is reported in Fig. |lj along with 
the role played by each vector component. The quark 



n 2 



"n 3 3 




FIG. 1. The Quark mixing torque vector 
mixing matrix can be written using the Cayley Hamilton 



3 



theorem @ (see Sect. HV|) as 

V = e rA \ T=1 = 1 + sinc(|ft|)A + i(sinc(^i)) 2 l 2 

| ft | = X v /l + y 2 X 2 +x 2 X i 

■ i \ sin(a) 

smc a = 

a 

(13) 

Moreover from Eq. (fT2"j) the action of the mixing ma- 
trix on the initial vector Q can be specified through the 
following Rodriguez rotation^ 

Q= cos(|Tf|)Q + sin(|Tt|)"^ x Q 

+(1 - cos(| ft \))(Jt ■ <3o) • It (14) 



n 

W\ 



The geometrical interpretation is less obvious if we in- 
clude the CP violating term. We assume Eq. (fTTj) to be 
still valid and with a slight abuse of the notation write 



ft = ft x + zft 2 



fli = (-yA 2 ,a;A 3 cos((5),-A) 
Tt 2 = (0,xA 3 sin(<5),0) 



(15) 



This assumption contains the bare essence of CP viola- 
tion from a geometrical point of view. The vector O 
splits into a real and imaginary part, as shown in Fig. 
@ where the second component of the torque vector is 
composed by two subcomponents: 

a) : the coupling vector ft 1,3 = (yA 2 ,0,— A) is the com- 

ponent of the vector in the 1-3 plane; 

b) : the CP violating sector is viewed as the pseudo vector 

(fJ 1>3) Im(0 2 ),Re(02)). 



b) 





FIG. 2. The real(a) and the imaginary part (b) of the torque 
vector Q. 

In terms of the Wolfenstein parameters the modulus of 
the Torque vector can be written as 



iTtl 



\ 



A 2 + (AX 2 ) 2 + 



(AX 3 ) 2 (16) 



|ft| = ^A2(l+tan 2 (0 A )) + A3 (p 



[l + tan 2 (5)} tan 2 (<^ 



(17) 



tan(0A) = AX, tan(5) = 



where 4>a and 8 are indicated in Figs. ©. 
The angle (f>A lies in the (1,3) sector and specifies the 
direction of the ft vector components in this plane. We 
visualize the geometric content of our problem as indi- 
cated in the second of Figs. ([2]), in which the complex 
vector component lying along the direction of the axis 2 
is split into an imaginary and a real part. 
In more rigorous mathematical terms we can illustrate 
the above procedure as it follows. We first note that 



A = A 1 +A 2 



(18) 



with 



A? 








A 3 \ 




1 









V- A 3 





: 





Ai 
A 2 
-A 2 



The matrices labelled with 2, 1 are not commuting 
each other, therefore we have at the first order in the 
Zassenhaus disentanglement formul43 [10] 



V = e A2+Al 



e A2 e Al e c 



C = -\[A 2 ,A l ] 
Where the (anti-hermitian) matrix C is given by 

C = 



(19) 




(20) 



Neglecting the matrix C, which is of the order o(A 4 ) 
find that the CKM matrix can be expressed as 



we 



V ^ e A2 e 



Ai a Ax 



The Zassenhaus formula writes e / 



(21) 



where 



the operators C m are given in terms of successive commutators, 
the first two being Ci = -5 [A, B], 62 = 5 [A, B] + [A, §]] 



4 



With 



V 2 



cos(|A 3 |) 




^|Sin(|A 3 |) 

1 



r^sin(|A 3 |) cosQAsI) 



(22) 



and the use of the Cayley Hamilton theorem allows the 
following (exact) form of the second exponential 



V 



( xX 3 e lS \ 


V -xX 3 e~' s 0/ 



/ \ 
yX 2 





( 


A 


0\ 






-A 





K 


V 








e Al = Vx = C t + C X A X + C 2 A\ 



1 



( C °) 

Ci 



( |Ai :2 | 2 

i 1-t -(1+t)/ 



|Ai, 2 | 2 
e *|Ai, 2 | 

fi -»|Al,2| 



V 2|A ll2 p / 
(23) 



T 



|A 1)2 |(l + »), |Ai, 2 | = ^/AfT 



The above formulae are a restatement of the tentative 
geometrical picture of Fig. ©. The nave disentangle- 
ment has reduced the CKM generation to the product of 
two matrices, V\ accounting for the mixing, induced by 
the vector S7 13 , and V 2 specifying a complex rotation, 
responsible for the CP violating contributions. 
The matrix (|2ip is an approximation of the exponen- 
tial form at the order o(A 4 ), but it is not equivalent to 
Wolfenstein matrix. The matrix (|21j). albeit an approxi- 
mation, since we have neglected higher order commuta- 
tors, is unitary at any order in the coupling parameter, 
while WVjy = 1 + o(A 4 ) (where Vw is the matrix (HJ)). 
We have stressed that the simple picture in terms of Eu- 
ler rotation is hampered by the presence of a complex 
term, the V matrix cannot be written in terms of the 
generators of rotations and indeed we find 



y _ e ~i(\R 1 +y\ 2 R 2 +x\ 3 T) 

e lS 
T = -i\ 

-e~ lS 



(24) 



The T matrix does not belong to SO (3) and the quark 
mixing matrix, written as the product of the exponential 
matrix correct up to the order o(A 4 ) is 



y _ e -tx\ 3 f e -iy\ 2 R 2 e -tX 3 Ri 



(A 4 ) 



C(xX 3 ) e lS S(xX 3 ) 

1 

-e- lS S{xX 3 ) C{xX 3 ) 

1 

C(yX 2 ) S(yX 2 ) 

-S(yX 2 ) C(yX 2 ) 

C((j)) = cos(0), S((p) = sm(4>) 



C(A) 5(A) 
-5(A) C(A) 
1 



and displays the largely well-known feature that the mix- 
ing angles are proportional to the Cabibbo coupling pa- 
rameter according to 

§i 3 cx A J = ./ 

V mi 



(26) 



^2,3 oc A = 



□ x ^ rn d 

1712 « A = . / — 
V TO., 



Furthermore, in full agreement with the particle data 
group paradigm, we get 

( xX 3 e' s \ 




V -xX 6 e~ l6 0/ 







1 



U s = e}\ with 5= 000 , R 3 = i\ 00 



s 



-10 



(27) 

It is also interesting to note that the T matrix can be 
written as 



T = -icos{5)R 3 +ism(5)S 3 

/0 1 
5 3 = [ 
10 0, 

The nave disentanglement (order o(A 4 )) 

y _ e -iR e ix\ 3 sin((5)S 3 

R= XR 1 + yX 2 R 2 + xX 3 cos(8)R 3 



(28) 



(29) 



5 



Corresponds to the product of two matrices, namely 
V S VrVi 



where 



V, = 



Vr 



C(xX 3 sm(S)) tS(x\ 3 sm(6)) 

1 

iS(xX 3 sm{5)) C(xX 3 sin(S)) 

C(AA 3 ?? ) -iS(AX 3 v ) \ 

10 
-zS^A 3 /?) C(ylA 3 r;) / 



= (30) 



and can be written as Eq. 1131 with 

\(-y\,x\ 2 cos(S), -1) 



(31) 



The matrix Vj mixes the first and third quark gen- 
eration mass eigenstates and is responsible for the CP 
violation. It is a pseudo rotation matrix and is generated 
by a matrix whose determinant is the Jarlskog invariant 
[111 ] . discussed in the forthcoming section. 
We have so far shown that the exponential parameter- 
ization implicitly contains Wolfenstein and Euler type 
forms, in the following sections we will dwell on its fur- 
ther advantages. 



IV. THE CAYLEY HAMILTON THEOREM 
AND THE QUARK MIXING MATRIX 

The exponential matrix ([T]) can be treated in different 
ways. 

We have already shown that the use of a Taylor expansion 
leads to a Wolfenstein form, which preserves the unitar- 
ity of V at the expansion order (the mixing matrix in Eq. 
© is unitary at the order o(A 4 )). 

The method of the exponential disentanglement can be 
used too and such a procedure allows an interesting ge- 
ometrical picture of the mixing dynamics and albeit an 
approximation in the Cabibbo coupling parameter, the 
mixing matrix written as in Eq. (|17[) preserves the uni- 
tarity at any order in A, as discussed more accurately in 
the concluding remarks. 

The matrix V can, however, be written in an exact form 
using the Cayley Hamilton theorem, by setting 



where 



V = C a l + CiA + C 2 A 2 (32) 



C + ejd + e 2 C 2 , with j = 1,2,3 (33) 



With Sj being the roots associated with the characteristic 
equation of the matrix A, namely 



in %- + iA = o 



(34) 



A = 2xyX 6 sm(S) = -2AtjX 6 
| it | = A v /l + y 2 A 2 +a; 2 A 4 = 



(35) 



= A,/l + (AA) 2 + (A 2 A4) (p-|r + 7? 



i A is the determinant of the matrix A. 

A little bit of algebra yields to define the C< (i = 0, 1, 2) 

coefficients as the product of two matrix 





( 


Ci 






V 



£l£3 



£l£2 



\ 



'(£2 + £3) 



-(£1 ^ 


-£3) - 


(£1 + 


1 




1 




e £l 




(£ 2 - 


£i)(£3 


-El) 








(£1- 


£2) (£3 


~£ 2 ) 




e £3 





£2 



(36) 



V (£1 - £3)(£2 - £3) / 

Eq. (|23p (along with Eqs. (|27p) is the most general form 
of the quark mixing matrix which can be derived from 
an exponential parameterization, it is exact but not easy 
to remember. 

Let us now give an idea of the orders of the numerical 
values characterizing the various quantities entering the 
above equations. The use of the previously quoted val- 
ues for the Wolfenstein parameters lead to the following 
evaluations for the solution of Eq. ([34]) 



ei = -0.231712 
e 2 = 0.00117*- 
e 2 = 0.23054i- 



(37) 



It is worth stressing that the matrix D provides the di- 
agonal forms of either V and A. It follows therefore that 
the two matrices have the same eigenvectors. They can 
be determined using A instead of V, because the proce- 
dure is significantly simpler. We find that the eigenvalues 
are in the form 



\3>- 



I 1 \ 

-Sj — xyX 5 e~ lS 

yX 3 — EjxX 3 e 



\ 



(38) 



/ 



It is worth mentioning the companion matrix asso- 
ciated with the characteristic equation (30) [12j, which 



6 



C(X) + 



V 



i(n- 



C(A) - (AA 2 ) 2 [i 



A /l 
2 3 \4 



IT) 



-ax s [(n - 1) - ^ (n - 1)] -s(aa 2 ) + ^1 (A 2 (n - ±) - i2(n - §)) 



AX 



[ir-£(ir-i; 



s(AX>) + Af (A 2 (n * - A) - i2(rr - |)) 



(39) 



(AA 2 ) 2 A 2 A 6 ^ 

1 I .1 I ^ 



-12(p 2 + j? 2 ) - 8iy + 12p- 2; n = p + i?7; IT = p - zt? 



writes 



/ -£i£2£3 \ 

1 -(£i£ 2 + £ 2 £3 +£l£3) 
1 -(£i+£ 2 +£ 3 ) 



(40) 



/ 



It is accordingly expressed in terms of three invariant^, 
namely 



£i£ 2 £3 = — 2ix\ sin(5) 
s 2 +s 2 +e 2 = \^\ 2 
£l + £2 + £3 = 



(41) 



the first of which is the Jarlskog invariant, a measure of 
the amount of CP violations, emerging in quite a natural 
way in the present analysis. 



V. CONCLUDING REMARKS 

We have shown that the exponential parameterization 
interpolates between Wolfestein and Euler like forms and 
could provide a useful and flexible tool of analysis. Its 
approximations in terms of the Cabibbo coupling can be 
either expressed as Taylor expansions or as unitarity pre- 
serving forms bassed on the Zassenhaus formula. 
The Taylor expansion does not meet too much aestheti- 
cal criteria, but it can usefully be exploited to get higher 
order approximations of Wolfenstein type parameteriza- 
tions an example is shown below, where we report the 
nave expansion of the exponential matrix up to the or- 
der o(A 7 ). 



4 A 3 X 3 matrix has three invariants given by its determinant, its 
trace and by the sum of the determinants of its minors 



We have reported the matrix (|39f (where C(A), S{\) de- 
note the expansion of cosine and sine up to the order 
o(A 7 )) for comparison purposes with other forms avail- 
able in literature. The accuracy of this last matrix is one 
part over 10 9 and can therefore considered exact for any 
expansion purposes. 

The extension of the CKM matrix to higher dimensions 
by the use of the exponential matrix method is not com- 
plicated. In the case of four quark generations, we define 
the matrix containing 2 CP violating phases, appended 
to the smallest coupling terms. We have furthermore as- 
sumed that the coupling strengths to the fourth family 
be of the order A 3+n , n = 1, 2, 3. 



A 



( A 

-A 

~e- lSl x\ 3 -yX 2 

\ -e~ lS2 zX 6 -pX 5 



e lSl xX 3 e lS2 zX 6 \ 



\ 3 e %&2 

yX 2 pX 5 

uA 4 

-uX 4 



(42) 



/ 



The relevant Wolfenstein like approximation of the mix- 
ing matrix is reported in (|43|) . 

Furtheremore the invariants (4 for a 4 x 4 matrix), 
obtained directly from (|42|) read 



J2 = 



J 3 



X 2 f(l,yX,xX 2 ) + X s f(u, P X,zX 2 ) 
where f(a, b, c) = a 2 + b 2 + c 2 

2iX & [xysin(Ji) + zpX 6 sin(<S 2 )J = 
xyX 6 (e lSl - e~ lSl ) + zpX 12 (e^ 2 - 



- tSr 



) 



J4 = u X + [2itj/zcos(<5 2 ) — 2pux cos(<5i)]A 



13 



[P 2 



- px 2 ye li - S2 ~ Sl) 



pxyze 



-t(«a-«i) 



xy 2 z]X 



Hi 



the first invariant, associated with the trace of A, is zero. 
It is evident that the J 2 and J3 invariants are just a gen- 
eralization of those reported in Eq. (|3"2"|) while the fourth 
is completely new being associated to the full determi- 
nant of the matrix. We have reported this example to 
show the flexibility of the method it is however evident 
that the detection of CP violating effects due to the new 
phase require an accuracy at least of the order A 6 . 
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V = 



C(A) + ^$ 

-s(x) - ^(n x - 1) 

-aa 3 [(n a - 1) - £ (Hi - !) 
-A 6 P (n 2 - 1) 



S(A) 



i(ni 



C(A)-(^A 2 ) 2 [| -£ (I-,,)] 
-S(AA 2 ) + ^ (A 2 ( ni - i) - 12(n a - D) 
—A 5 (_ P + Au|) 



ylA' 



[n; 



-n^pA 6 \ 



5(aa 2 ) + ^- (A 2 (n^ - A) _ i2(n; - §)) a 5 ( p + Au|) 



+ 



2 1 4! 

-kA 4 



uA 4 
1 



X 1 

$1 = pi + ir)i with pi = - cos(Ji) + -, T]i = 
2/ 2 



(43) 



X z \ z 

--sin(5i); $2 = pi + ir)2 with p 2 = - cos^) + -, 7/2 = — sin(<5 2 ) 



Pj - Mfe J 



1,2 



Before concluding the paper we will address the prob- 
lems associated with the exponential forms of the neu- 
trino mixing matrix, which have also been discussed in 
fl3} , where the lepton-quark complementarity [lij has 
been reformulated by noting that the relevant rotation 
occur around axes forming an angle of 45°. The present 
experimental data seem to favor the tribimaximal (TBM) 
form [l5j therefore the neutrino mixing matrix reads 



tries of the B matrix and those of the TBM form 

1 \ 



B = a 







V 



V2+1 



V2 + 1 



u 



[2 


1 


\ 




V 3 






1 


1 


1 




"71 


71 


71 


(44) 


1 


1 


1 




71 


"71 


71 





If we assume that also this form is generated by an ex- 
ponential matrix (with all real entries) according to 



/ 



U = e u : B 



a (3\ 
7 

-7 

/ 



(45) 



We obtain the following correspondence between the en- 



a 



V2V2 + 2V6 + 9 



sm 



(46) 

The values of the entries of the TBM matrix do not allow 
the interpretation of the neutrino mixing matrix as an ex- 
pansion around the unit, notwithstanding it is possible 
to get a better agreement with experimental by making 
an appropriate expansions around the matrix B and then 
around the TBM, as it will be shown in a dedicated pa- 
per. 

In this paper we have provided an extensive account of 
the possibilities offered by the exponential form of the 
CKM matrix, which looks like a prototype from which 
all the other forms can be derived, we hope that our sug- 
gestions provide a useful tool in the relevant applications. 
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